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SUMMARY 

An  analytical  procedure  has  been  developed  to  predict  the  circumferential 
pressure  profile  at  the  exit  of  a  compressor  stage  subjected  to  a  spatially 
steady  inlet  distortion.  Expressions  have  been  developed  relating  the  pressure 
ratio  and  weight  flow  characteristics  of  an  axial- flow  compressor  stage  to  the 
normal  force  and  drag  characteristics  of  an  isolated  airfoil.  These  steady- 
state  interrelations  between  rotor  and  isolated  airfoil  were  used  to  apply 
available  unsteady  data  for  isolated  airfoils  in  the  study  of  distorted 
inflow  effects.  Both  unsteady  and  quasi- steady  predictions  were  made  and 
were  compared  with  available  experimental  results  from  low- speed  compressor 
tests.  It  was  found  that,  in  general,  the  unsteady  predictions  are  in  better 
agreement  with  experiment  than  the  quasi- steady  predictions,  but  there  are 
regions  where  both  theories  are  in  significant  disagreement  with  the  data. 

The  results  of  the  work  reported  herein  have  led  to  an  improved  analytical 
model  which  is  briefly  reviewed.  Further  analytical  work,  based  on  this 
improved  model,  will  be  the  subject  of  future  reports. 


INTRODUCTION 


The  performance  of  an  axial-flow  compressor  can  be  seriously  compromised 
by  the  presence  of  a  nonuniform  circumferential  velocity  distribution  in  the 
inlet  of  the  compressor.  Such  an  inlet  velocity  distortion  could  be  caused 
by  a  poorly  designed  engine  inlet,  flow  separation  from  some  portion  of  the 
inlet  cowl  arising  from  a  locally  large  incidence  angle,  or  even  by  scavenging 
of  a  thick  boundary  layer  from  some  adjacent  aerodynamic  surface.  Whatever 
the  cause,  the  primary  effect  will  generally  be  the  same;  viz.,  if  there 
exists,  circumferentially,  some  region  of  reduced  axial  velocity,  a  compressor 
blade  travelling  at  a  constant  rotational  speed  will  experience  a  periodic 
change  in  incidence  angle,  as  shown  in  Fig.  1.  The  aerodynamic  performance 
of  a  blade  in  a  distorted  flow  will  depart  from  that  for  some  nominal  operat¬ 
ing  condition;  the  degree  of  this  departure  will  depend  largely  on  such 
parameters  as  the  number  of  distortions  per  revolution,  the  c  „’cumferential 
extent  of  each  distortion,  the  magnitude  of  each  distortion,  and  the  rate  at 
which  incidence  angle  changes  due  to  the  distortion  taJc-j  place.  A  change  in 
compressor  blade  performance  from  some  nominal  value  will  affect  the  overall 
performance  of  the  compressor  stage.  In  many  cases  a  distortion  is  severe 
enough  to  affect  the  performance  of  the  entire  engine  to  the  extent  that  the 
engine  surge  margin  will  be  compromised. 

Current  techniques  for  predicting  engine  stability  with  a  distorted 
inflow  treat  the  problem  as  quasi- steady  and  thus  employ  steady- state  stage 


characteristics.  These  techniques  assume  that  the  compressor  blading  responds 
instantaneously  to  distortion- Induced  changes  in  the  incidence  angle.  The 
response  of  a  compressor  blade  to  a  rapid  change  in  inflow  is,  in  fact,  a 
nonsteady  phenomenon  that  is  far  from  instantaneous.  This  is  indicated  by 
the  results  shown  in  Fig.  2  which  were  taken  from  Ref.  1.  In  this  study, 
circumferential  distortions  having  a  total  circumferential  extent  of  180  deg 
were  configured  in  three  possible  ways:  (a)  a  single  180  deg  distortion, 

(b)  two  radially  opposed  90  deg  distortions,  or  (c)  four  equally  spaced 
45  deg  distortions.  This  figure  shows  that  the  surge  pressure  ratio  (plotted 
vertically)  is  strongly  affected  by  the  change  in  circumferential  extent  of 
each  portion  of  the  distortion,  even  though  the  total  circumferential  extent 
of  the  distortion  remained  invariant.  (This  observation  is  further  sub¬ 
stantiated  by  the  results  shown  in  Fig.  2  of  Ref.  2  in  which  the  total 
circumferential  extent  of  the  distortion  was  90  deg. ) 

The  sensitivity  of  compressor  characteristics  noted  above  to  the  time 
rate  of  change  of  the  distortion  pattern  (as  observed  in  the  relative  system 
of  the  blades)  was  anticipated  by  the  analysis  of  unsteady  aerodynamic  blade 
response  to  variations  in  angle  of  attack  described  in  Ref.  3.  In  this  work 
use  was  made  of  a  body  of  unsteady  aerodynamic  data  for  an  isolated  airfoil 
oscillating  both  in  potential  flow  and  well  into  the  stall  regime.  By 
applying  these  unsteady  aerodynamic  data  to  a  hypothetical  distortion,  the 
analysis  examined  the  response  of  a  compressor  blade  to  several  distortion 


variables,  including  time  rate  of  the  distortion  pattern  and  degree  of 
penetration  into  the  stall  regime.  This  preliminary  analysis  predicted  the 
expected  enhancement  in  lift  and,  within  its  limited  framework,  implied  the 
resulting  pressure  rise  that  can  result  from  certain  distortion  patterns,  as 
«ell  as  the  flow  breakdown  from  other  patterns. 

The  current  investigation  is  directed  toward  assessing  the  exrtent  to 
which  the  characteristics  of  compressors  operating  in  distorted  flow  are 
influenced  by  unsteady  aerodynamic  effects.  In  this  analysis,  refinements 
have  been  made  to  the  techniques  of  Ref.  3  to  permit  a  more  complete  prediction 
of  the  characteristics  of  a  compressor  in  an  arbitrarily  distorted  inlet  flow. 
In  addition  to  refining  the  unsteady  lift  prediction  developed  in  Ref.  3>  a 
method  for  converting  these  results  to  stage  pressure  rise  has  also  been 
developed,  and  results  have  been  compared  with  available  experimental  data. 

This  report  covers  the  results  obtained  during  the  first  year  of  this 
contract.  This  includes  the  derivation  of  the  basic  fluid  dynamic  equations 
and  the  comparison  between  the  predicted  and  measured  unsteady  compressor 
performance.  Subsequent  reports  will  cover  the  continuation  of  this  analysis 
including  an  examination  of  the  assumptions  made  in  this  first  study  and 
refinements  made  to  the  theory. 
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The  problem  under  consideration  in  this  report  involves  the  nonuniform 
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flow  of  fluid  through  a  compressor  blade  row.  In  its  most  general  form  this 
problem  is  virtually  impossible  to  solve,  and  it  is  necessary  to  reduce  its 
scope  considerably  to  make  it  tractable.  One  essential  characteristic  of 
the  problem  is  that  a  compressor  blade  row  is  moving  at  constant  rotational 
speed  through  a  fluid  medium  having  a  variable  vector  velocity.  As  described 
in  the  Introduction,  this  produces  a  time-varying  angle  of  attack  on  the 
individual,  compressor  blades,  and  it  is  this  essential  nonsteady  character¬ 
istic  that  will  be  studied  in  the  present  report,  although  other  elements 
of  the  problem  will  be  simplified  for  clarity. 

In  the  analysis  which  follows,  a  number  of  fundamental  assumptions  have 
been  made  to  simplify  this  first  effort  at  finding  a  solution.  These  assump¬ 
tions  are  listed  below. 

1.  The  flow  is  inviscid  and  incompressible. 

2.  The  flow  through  each  blade  passage  is  quasi- steady. 

3.  The  axial  velocity  at  any  circumferential  station  is  constant  through 
the  stage.  ' 
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5.  The  flow  through  adjacent  blade  passages  is  periodic. 

6.  Three-dimensional  radial  effects  are  negligible  although  radial 
equilibrium  corrections  are  applied  to  the  measured  static  pressure. 

7.  Unsteady  effects  manifest  themselves  at  the  exit  plane  through  the 
unsteady  aerodynarJ  c  response  of  the  blading  only. 

8.  Jnsteady  stalling  on  a  cascade  blade  can  be  modeled  from  the 
unsteady  stalling  response  of  an  isolated  blade. 

9.  Unsteady  drag  can  be  computed  as  a  departure  from  the  steady  drag, 
proportional  to  the  unsteady  lift  departure  from  steady  lift. 

Quasi- Steady  Flow  Through  a  Blade  Row 

The  blade  row  to  be  analyzed  here  is  shown  in  Fig.  3  together  with  some 
of  the  nomenclature  to  be  used.  Station  1  is  located  at  the  blade  row 
entrance  plane  and  station  2  is  located  at  the  exit  plane.  The  tangential 
blade  gap  is  r  and  the  blade  chord  is  c  =  2b.  The  blade  is  oriented  at  a 
fixed  chord  angle,  a^,  relative  to  the  tangential  direction.  Axial  quantities 
are  denoted  by  the  subscript  A  and  tangential  quantities  by  the  subscript  Q,. 

The  two  right-hand  velocity  triangles  are  drawn  for  the  absolute  reference 
frame,  stationary  in  space,  and  consist  of  the  absolute  velocity,  V  ,  the 
wheel  speed,  U,  and  the  relative  velocity,  V.  These  quantities  are  sub¬ 
scripted  1  and  2  for  entrance  and  exit,  respectively.  The  left-hand  velocity 
triangles  are  included  for  convenience  and  denote  the  quantities  of  major 
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interest  in  the  relative  frame  of  reference  which  are  the  relative  velocity, 


V,  the  axial  velocity,  VA>  and  the  tangential  swirl  component,  V<^.  The  total 
force  on  each  blade  is  denoted  by  F  and  is  resolved  into  axial  and  tangential 
components,  FA  and  F^,  or  into  lift  and  drag  components,  L  and  D.  For  this 
analysis  the  lift  is  assumed  to  act  normal  to  the  inlet  relative  velocity 
vector,  V-j_. 

In  Fig.  4  the  inlet  and  exit  velocity  triangles  are  shown  in  more  detail. 
The  absolute  components  are  shown  as  solid  vectors  and  the  axial  and  tangen¬ 
tial  swirl  components  are  shewn  as  clashed  vectors.  In  this  diagram  all 
angles  are  measured  relative  to  the  tangential  direction.  The  angle  of  the 
relative  velocity  is  denoted  by  (3,  and  that  of  the  absolute  velocity  is 
denoted  by  a.  In  addition,  the  yaw  angle,  ay>  of  the  inlet  absolute  velocity 
is  shown  on  the  inlet  triangle.  (This  angle  is  the  supplement  toa^.) 

For  a  given  inlet  distortion  the  circumferential  variations  of  the  inlet 
total  and  static  pressures,  P01  and  p^,  and  of  the  inlet  yaw  angle,  ay,  will 
be  given  quantities,  usually  obtained  from  an  experiment.  In  keeping  with 
the  assumptions  listed  in  a  previous  section,  the  steady  form  of  Bernoulli's 
equation  will  be  used  to  confute  the  instantaneous  velocity  magnitude  into 
the  blade  row  at  any  circumferential  position,  0  , 

P,<*>  +  2  />V*,(0)=  Po,(0)  (1) 


dagaaaaiiaaaaffiffeaaia  -  r—  • aaaifiiStiato aflgasaaaaa£ -  -  -  - 


Note  that  the  density  is  unsubscripted  in  accordance  with  the  assumption  that 
the  flow  is  incompressible.  This  equation  may  be  solved  now  for  the  absolute 


inlet  velocity. 


p„,(e>-  p,<9). 


The  aerodynamic  force  response  of  the  rotor  blades  will  depend  on  the 
circumferential  variation  in  the  inlet  air  angle,  0^.  From  the  inlet  velo¬ 
city  triangle  of  Fig.  4  it  is  easily  shown  that 


and  hence 


VA  (0)  =  V6|(0)  sin  Oy (0) 


Vn(0):  U  +  VO({0)  cos  ety  (0) 


VA  (0)  V  (0)  sin  a  (0) 

/3'(0)  (0 )  U  +  Vo,(0)cosay  (0) 


It  is  more  convenient  to  calculate  blade  aerodynamic  response  from  the  rotor 
blade  angle  of  attack,  a^,  which  is  shown  in  the  upper  half  of  Fig.  5  to  be 


The  procedure  to  be  followed  in  this  analysis  is  as  follows.  Prom  the 
circumferential  variation  of  the  inlet  pressures  and  yaw  angle  we  can  cal¬ 
culate  V01  from  Eq.  (2),  from  Eq.  (5)  and  aR  from  Eq.  (6).  If  the  com¬ 
pressor  wheel  rotates  at  a  constant  rate  SI  =  d9/dt  then  0  =  fit  and  these 
quantities  become  time  histories.  From  the  time  history  of  a.  it  iB 

ft 

possible  to  enter  a  table  of  unsteady  aerodynamic  lift  and  drag  data  appropriate 
to  the  rotor  blade  and  obtain  unsteady  lift  and  drag  time  histories  versus  9  . 
Details  of  this  procedure  will  be  given  in  a  subsequent  section.  Next  the 
lift  and  drag  forces  can  be  resolved  into  axial  and  tangential  force  components 
(see  bottom  of  Fig.  5)  from  the  equations 


Fa  (9)--  L(0)co»  /3,(0)  -  D(0)  sin 
^(0)  =  L(e)S!n^,(a)  +  D(fi)  cos /3,(0) 


(7) 


Further  details  of  the  quasi- steady  through- flow  analysis  are  given  in 
Appendix  I.  Use  is  made  of  the  continuity  and  momentum  equations  to  obtain 
formulas  for  the  static  and  total  circumferential  pressure  distributions  at 
the  rotor  exit  plane,  which  are  repeated  below  for  convenience. 


/>V2  (0) 

pJ0)  =  p.(0)+  -T-. -  C  (Q) 

2  i  2r/c 


(8) 
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Use  of  Unsteady  Aerodynamic  Data 


In  attempting  a  solution  of  the  unsteady  inlet  distortion  problem,  the 
unsteady  lift  and  drag  data  which  should  be  used  should  ideally  be  obtained 
for  blades  in  cascade  operating  into  the  stall  regime.  Unfortunately  such 
data  are  net  available  for  cascades,  and,  moreover,  there  is  no  known  body 
of  unsteady  drag  data  obtained  on  any  lifting  surface.  An  attempt  has  been 
made  in  the  present  analysis  to  overcome  these  deficiencies  in  the  data  by 
using  available  unsteady  lift  data  on  an  isolated  airfoil  and  by  hypothesizing 
a  reasonable  unsteady  drag  model.  Suitable  transformations  have  been  made  to 
convert  these  isolated  airfoil  forces  into  their  equivalent  cascade  counter- 
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total  lack  of  unsteady  cascade  data,  and  the  analysis  is  so  structured  that 
if  cascade  data  were  to  become  available  it  would  be  a  simple  matter  to  insert 
it  in  place  of  the  isolated  airfoil  data. 

A  number  of  steps  are  required  to  introduce  the  available  unsteady  lift 
data  for  the  isolated  airfoil  into  the  pressure  equation  to  represent  the 
unsteady  lift  on  a  rotor  blade.  These  steps  are  summarized  briefly  in  this 
paragraph  to  indicate  how  all  of  the  parts  fit  together,  and  the  details  of 
each  step  are  discussed  in  Appendix  II.  First  the  equations  must  be  derived 
that  relate  steady- state  cascade  blade  lift  to  compressor  blade  performance. 
This  is  done  by  manipulating  Eqs.  (7),  (8),  and  (9)  together  with  other  rele¬ 
vant  equations.  Next  a  scaling  law  must  be  determined  that  will  relate  the 
lift  on  a  cascaded  airfoil  to  the  lift  on  the  isolated  airfoil  (see  Appendix 
III) .  The  unsteady  isolated  airfoil  lift  must  then  be  introduced  and  scaled 
to  yield  an  incremental  change  in  the  unsteady  rot or  blade  lift  which  will 
subsequently  be  used  to  calculate  the  unsteady  rotor  pressure  response.  A 
similar  procedure  for  scaling  unsteady  drag  and  details  of  the  unsteady  drag 
hypothesis  are  in  Appendix  IV. 

A  complete  summary  of  all  of  these  steps,  including  the  sequence  of 
computer  operations,  is  in  Appendix  V. 
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APPLICATION  TO  A  COMPRESSOR  ROTOR 

Description  of  UARL  Three- Stage  Compressor  Rig 

The  brief  description  which  follows  of  the  three- stage  compressor  rig  is 
included  in  the  present  report  to  assist  the  reader  to  understand  the  experi¬ 
mental  data  discussed  later.  It  does  not  purport  to  be  a  detailed  description 
of  this  facility,  nor  does  it  attempt  to  provide  extensive  coverage  of  the 
experimental  results  obtained. 

The  three-stage  rig  is  a  low- speed,  low-pressure  facility  consisting  of 
an  inlet  bellmouth,  inlet  guide  vanes  ( IGV) ,  and  three  conventional  rotor- 
stator  stages.  The  rotating  parts  were  driven  by  a  500  hp  dc  electric  motor. 
The  airflow  passage  has  an  outer  diameter  of  24  inches,  and  the  hub- to- tip 
ratio  for  the  tests  described  herein  was  0.7.  A  centerline  section  drawing 
is  presented  in  Fig.  6,  approximately  to  scale.  The  rotors  are  denoted  by 
Rl,  R2,  R3,  and  the  stators  by  SI,  S2,  S3.  The  instrumentation  stations  are 
numbered  above  the  drawing,  and  the  pressure  instrumental ion  used  at  each 
station  is  shown  in  the  sketches  at  the  bottom  of  the  figure.  Total  pressure 
measurement  locations  are  denoted  by  the  circled  points  and  static  pressure 
measurement  locations  by  the  diamond  points.  In  addition,  yaw  probes  (not 
shown)  were  used  to  make  flow  angle  measurements  at  the  entrance  to  each 
rotor.  The  distortion  screen  was  located  in  a  constant  area  region  of  the 


annulus  ahead  of  the  IGV  and  consisted  of  an  annular  segment  of  thin  perforated 
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plate  (approximately  50  percent  porosity).  The  segment  angle  of  most  concern 
in  the  present  study  was  180  deg,  although  tests  were  also  run  with  distortion 
screens  of  90  deg  and  135  deg  angular  extent. 

The  normal  test  procedure  was  to  first  obtain  a  baseline  characteristic 
with  no  distortion  by  varying  the  downstream  throttle  and  making  pressure 
measurements  at  a  constant  speed  condition.  This  is  shown  in  Fig.  7  (and  is 
described  in  Appendix  II)  for  the  first  stage.  A  distortion  screen  was  then 
inserted  and  measurements  were  made  at  the  same  rotor  speed  and  at  some  fixed 
downstream  throttle  setting,  representative  of  a  given  undistorted  point  on 
the  characteristic  curve.  Because  there  were  only  eight  circumferential 
measuring  stations  for  each  axial  station,  it  was  necessary  to  repeat  this 
test  a  number  of  times  with  the  distortion  screen  incremented  circumferentially 
each  time.  Typically,  a  total  of  10  or  more  test  points  might  be  taken  to 
obtain  the  desired  circumferential  pressure  variations. 

Detailed  Study  of  a  Single  Distorted  Flow  Condition 

In  the  performance  of  this  work  a  number  of  distorted  flow  conditions 
were  studied,  and  many  variants  to  the  analysis  were  introduced,  both  to 
improve  the  correlation  between  analysis  and  test  and  to  provide  a  deeper 
insight  of  the  mecnanisms  involved.  In  this  report  we  shall  consider  one 
of  these  distorted  flow  conditions  in  considerable  depth.  In  addition,  a 
mere  general  discussion  will  follow  in  a  later  section  of  other  cases  con¬ 
sidered  in  this  study,  and  of  the  sensitivity  of  the  results  to  variations  in 
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the  inlet  parameters.  In  all  cases  the  analysis  was  confined  to  the  first 
stage  only. 

In  this  detailed  discussion  we  shall  he  concerned  with  a  180  deg  distor¬ 
tion  of  the  inlet  flow  with  the  compressor  operating  at  the  peak  pressure 
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point.  If  no  distortion  had  been  present  this  point  would  have  been  near 
the  maximum  point  on  the  rotor  4>  -  0  curve  in  Fig.  7«  The  measured  circum¬ 
ferential  mid-annulus  total  pressure  profiles  at  the  rotor  inlet  and  exit 
planes  are  shown  in  Fig.  8,  the  outer  diameter  static  pressure  profiles  at 
the  rotor  inlet,  rotor  exit,  and  stator  exit  planes  are  shown  in  Fig.  9> 
and  the  yaw  angle  profile  at  the  rotor  inlet  plane  is  shown  in  Fig.  10.  The 
total  pressure  was  measured  in  a  stationary  frame  of  reference  so  there  is 
no  significant  change  in  total  pressure  from  rotor  exit  to  stator  exit.  The 
pressure  profiles  are  plotted  relative  to  the  local  atmospheric  total  pressure 
which  in  this  case  is  Pref  =  4o4.33  in,  HgO.  The  position  of  the  screen,  which 
extended  from  0  =  270  deg  to  0=  90  deg,  is  indicated  by  the  shaded  region 
along  the  abscissa  of  these  and  subsequent  figures. 

The  procedures  outlined  above  and  in  Appendix  V  were  followed,  and  the 
pressure  difference  P01(8)  -  p^e),  together  with  the  yaw  angle  variation, 
were  used  to  calculate  the  dimensionless  stall  angle  parameter,  a  ,  shown 
in  Fig.  11  (where  o  =-l  is  above  the  stall  angle  and  a<l  is  below  the  stall 
angle).  It  is  seen  that  the  rotor  blade  operates  above  stall  over  most  of 
the  circumference,  particularly  in  the  regions  behind  the  screens.  This 
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distribution  was  numerically  differentiated  to  obtain  the  dimensionless 
angular  velocity  and  angular  acceleration  parameters,  A  and  B,  as  described 
in  Appendix  III,  and  unsteady  lift  and  drag  data  on  the  rotor  blade  were 
obtained.  These  are  shown  in  Figs.  12  and  13,  respectively,  where  the  un- 
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steady  result  (solid  line)  is  cougar ed  with  the  quasi- steady  result  (dashed 
line) .  The  latter  was  obtained  by  setting  A  =  B  =  0  and  performing  the 
analysis. 

It  can  be  seen  that  there  are  large  differences  between  the  quasi-steady 
and  the  unsteady  results  over  much  of  the  range  of  operation  in  both  Figs.  12 
and  13.  The  only  region  where  the  quasi- steady  and  the  unsteady  results  agree 
is  the  circumferential  range  2OO<0<24O  deg.  In  that  region  a  has  a  con¬ 
sistently  small  gradient,  as  shewn  in  Fig.  11.  Very  large  fluctuations  occur 
near  the  edges  of  the  screen,  particularly  for  the  lift.  This  is  most  likely 
artificially  caused  by  the  magnification  of  error  that  occurs  when  a  function 
that  is  not  perfectly  smooth  is  numerically  differentiated.  In  the  analysis 


and  discussion  that  follows  these  effects  are  discounted,  and  we  feel  it  will 
be  possible  to  modify  the  analytical  model  in  the  future  to  overcome  these 
wide  excursions. 

Let  us  consider,  in  greater  detail,  the  large  differences  between  the 
quasi-steady  and  the  unsteady  results.  As  suggested  above,  we  will  exclude 
from  the  present  discussion  uke  regions  of  very  rapid  change  near  the  screen 
edges.  In  the  region  behind  the  screen,  over  the  range  -40<9'=40  deg  (or. 
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using  the  0  values  in  Figs.  11,  12,  and  13,  320  deg<0  and  0^40  deg),  a  is 
a  steadily  increasing  function,  and  in  this  region  the  unsteady  lift  is 
greater  than  the  quasi- steady  lift,  and  the  unsteady  drag  is  smaller  than 
the  quasi- steady  drag.  Conversely,  in  the  region  of  no  screen,  over  the 
range  115  <  0^200  deg,  a  is  a  steadily  decreasing  function,  and  in  this 
region  the  unsteady  lift  is  smaller  than  the  quasi- steady  lift,  and  the 
unsteady  drag  is  greater  than  the  quasi- steady  drag.  The  enhancement  of 
the  unsteady  lift  for  increasing  angle  of  attack  and  the  deterioration  of 
lift  for  decreasing  angle  of  attack  are  to  be  expected  in  this  high  angle 
of  attack  regime  (Ref.  3).  The  noted  changes  in  unsteady  drag  appear  to  be 
reasonable  when  related  to  the  lift  behavior.  These  differences  have  a 
great  deal  of  bearing  on  the  prediction  of  the  unsteady  exit  pressures  for 
the  rotor. 

The  circumferential  rotor  lift  and  drag  variations  were  next  trans¬ 
formed  to  axial  and  tangential  force  variations,  which  were  then  inserted 
into  Eqs.  (8)  and  (9)  to  predict  the  circumferential  static  and  total  pres¬ 
sure  profiles  at  the  rotor  exit.  This  was  done  for  both  the  quasi- steady 
and  the  unsteady  conditions  and  the  results  are  shown  in  Figs.  14  and  15  for 
the  total  and  static  pressures,  respectively.  The  parameters  used  in  the 
analysis  leading  up  to  Eqs.  (8)  and  (9)  were  obtained  at  mid-annulus,  and 
hence  the  predicted  pressure  profiles  in  Figs.  14  and  15  also  represent  mid¬ 
annulus  values.  Tne  exit  total  pressures  were  measured  at  mid- annulus,  and 
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are  shown  in  Fig.  14  as  circled  points.  However,  the  measured  static  data 
were  obtained  from  sidewall  taps  on  the  outer  diameter  and  had  to  be  corrected 
by  a  radial  equilibrium  criterion  to  yield  mid- annulus  results.  This  was 
done  using  the  formula  derived  in  Appendix  VI  and  the  results  are  plotted 
as  the  circled  points  in  Fig.  15. 

First  consider  the  total,  pressure  prediction  in  Fig.  14.  The  unsteady 
;heory  is  in  better  agreement  with  the  experimental  data  than  the  quasi¬ 
steady  theory  over  most  of  the  range  0  < 9-^270  deg.  The  sole  exception  here 
is  the  exit  from  behind  the  screen  in  the  vicinity  of  9  =  90  deg.  At  the 


entrance  to  the  screen,  over  the  range  270  <  9  <  360  deg,  both  theories  are 
in  substantial  disagreement  with  experiment.  The  static  pressure  comparison 
is  shown  in  Fig.  15  and  again  it  is  seen  that  over  much  of  the  circumference 
(with  the  exception  of  the  two  screen  edges  and  a  small  region  surrounding 
9=0  deg)  the  unsteady  theory  is  consistently  in  better  agreement  with  the 
measured  results  than  is  the  quasi-steady  prediction.  Furthermore,  over 
much  of  the  range  excellent  agreement  is  achieved  using  the  unsteady  theory. 
However,  this  agreement  tends  to  be  masked  by  th*.  large  changes  in  both  the 
unsteady  and  quasi-steady  pressure  profiles  due  to  the  circumferential,  changes 
in  the  inlet  pressures. 

Some  of  these  masking  effects  can  be  eliminated  by  plotting  the  rotor 


static  pressure  coefficient,  0  R)  as  a  function  of  9.  This  was  computed 
from  Eq.  (40)  in  Appendix  II  using  the  measured,  the  quasi-steady,  and  the 
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unsteady  circumferential  inlet  and  exit  static  pressure  profiles,  p1(9)  and 

P2 ( 8 ) •  The  results  are  shown  in  Fig.  16.  It  is  seen  here  that  the  unsteady 

prediction  agrees  reasonably  well  with  the  measured  ^  over  most  of  the 

K 

circumference  with  improved  agreement  over  the  range  110<9<220  deg.  In 
contrast,  the  quasi- steady  prediction,  shown  by  the  dashed  line,  has  the 
correct  order  of  magnitude  but  fails  to  predict  the  details  of  the  measured 
static  pressure  variations.  (Disagreement  between  the  unsteady  prediction 
and  measured  data  occurs  at  the  screen  edges,  caused  by  the  sudden  changes 
in  the  numerically  determined  derivatives,  as  explained  earlier  in  connection 
with  Fig.  12.) 

Additional  Distortion  Profiles 

Other  cases  were  analyzed  in  the  course  of  this  investigation.  Two 
additional  distortion  profiles  of  180  deg  extent  were  considered  which  are 
closely  related  to  the  case  discussed  above.  The  first  of  these  was  for  a 
higher  mass  flow  and  lower  pressure  ratio  than  the  peak  pressure  condition 
just  examined  and  is  denoted  as  the  midflow  condition.  The  nominal  undis- 
torted  steady- state  operating  condition  for  this  case  is  to  the  right  of  the 
peak  value  of  in  Fig.  7.  The  second  of  these  was  for  a  highly  loaded 
condition  and  is  denoted  as  operation  near  the  surge  point.  The  nominal  undis¬ 
torted  steady- state  operating  condition  for  this  case  is  to  the  left  of  the 

peak  value  of  £  in  Fig.  7. 

R 


The  measured  inlet  and  exit  mid- annulus  total  pressure  profiles,  outer 
diameter  static  pressure  profiles,  and  inlet  yaw  angle  profile  for  the  mid¬ 
flow  condition  are  shown- in  Figs.  17,  18,  and  19,  respectively.  This  is 
followed  hy  a  comparison  between  the  predicted  and  measured  static  pressure 
coefficient  for  the  same  case  in  Fig.  20.  A  similar  set  of  figures  for 
operation  near  the  surge  point  is  found  in  Figs.  21  through  24.  In  Fig.  20 
for  the  midflow  condition  we  have  fair  to  good  agreement  between  unsteady 
theory  and  experiment  behind  the  screen  (from  0  =  290  deg  to  0s?8o  deg)  and 
good  agreement  upon  emergence  from  the  screen  ( over  the  range  140 <fl< 240  deg) . 
The  same  statements  may  be  made  for  operation  near  the  surge  point,  shown  in 
Fig.  24,  where  the  best  agreement  is  found  over  the  ranges  20  <  0«=  90  deg  and 
160  c 0  < 350  deg.  In  ail  of  these  cases  exceptions  must  be  made  for  regions 
of  excessive  scatter  occasioned  by  the  numerical  differentiation  process. 

Additional  cases  were  examined  for  a  135  deg  distortion  at  peak  pressure 
and  for  a  90  deg  distortion  near  stall.  Once  again  the  measured  total  pres¬ 
sures,  static  pressures,  yaw  angles  and  unsteady  pressure  predictions, 
plotted  as  a  function  of  9,  will  be  found  in  Figs.  25  through  28  for  the 
135  deg  distortion  and  in  Figs,  29  through  32  for  the  90  deg  distortion. 

In  Fig.  28,  for  the  peak  pressure  condition,  both  the  unsteady  prediction 
and  the  quasi- steady  prediction  are  in  poor  agreement  over  most  of  the 
circumference.  In  Fig.  32,  for  operation  near  stall,  the  agreement  between 
unsteady  prediction  and  measured  results  is  fair  to  good  in  the  region  of 
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blade  emergence  from  the  distortion  screen  with  the  exception  of  seme  sharp 
discrepancies  at  the  screen  edge.  This  generally  good  agreement  persists 
over  the  range  90  <9<:230  deg.  Beyond  this  point  the  agreement  is  poor. 

A  critical  examination  of  all  five  cases  (Figs.  16,  20,  24,  28,  32)  shows 
that  the  best  agreement  between  unsteady  prediction  and  measured  response 
occurs  in  the  vicinity  of  blade  emergence  from  the  distortion  screen,  and 
poor  agreement  exists  in  the  vicinity  of  blade  entry  behind  the  distortion 
screen.  From  the  evidence  available  it  appears  that  the  unsteady  model  is 
most  effective  in  characterizing  the  recovery  from  stall,  with  good  agreement 
between  the  prediction  and  the  measured  data.  Further,  it  appears  that  the 
quiet  flow  Just  prior  to  entry  into  the  distortion,  and  the  disturbed  flow 
during  and  after  distortion  entry  is  more  or  less  poorly  modeled  by  the 
unsteady  analysis.  However,  before  too  many  generalizations  are  made  it 
must  be  realized  that  the  model  as  it  is  described  herein  is  relatively 
primitive  and  is  therefore  limited  in  its  capabilities.  Work  currently 
in  progress  has  suggested  that  it  may  be  possible  to  modify  this  model  to 
obtain  a  significant,  improvement  in  the  agreement  between  theory  and  experi¬ 
ment.  This  work  will  be  the  subject  of  a  future  report  under  the  current 
subcontract. 
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Stage  Pressure  Rise 

The  results  discussed  above  have  been  restricted  to  the  rotor  exit  plane, 
or  station  2  in  Figs.  3  and  6.  However,  it  is  desirable  to  be  able  to  pre¬ 
dict  (or  measure)  +he  stage  static  pressure  rise  which  requires  that  the 
static  pressure  through  the  stator  also  be  obtained.  (The  total  pressure  is 
referred  to  a  stationary  reference  frame  and  does  not  change  between  the 
rotors.)  Calculation  of  static  pressure  through  the  stator  is  a  relatively 
straightforward  procedure  if  the  stator  steady- state  characteristics  are 
known  in  advance,  or  can  be  accurately  predicted. 

In  the  present  case  the  stator  characteristics  were  available  from 

prior  measurement,  and  are  shown  in  Fig.  33  as  ^  versus  stator  inlet  angle 

S 

a 2  (see  Figs.  3  and  4).  The  procedure  to  be  followed  here  is  outlined 
briefly  below  and  is  described  in  detail  in  Appendixes  I  and  V.  The  rotor 
exit  velocity  triangle  is  used  to  compute  the  circumferential  variation  in 
stator  inlet  angle,  a2.  At  each  9  station  Fig.  33  is  entered  at  the  appro¬ 
priate  value  of  a  and  the  stator  pressure  coefficient,  0  ,  is  obtained. 
Equation  (40)  in  Appendix  II  is  then  rewritten  for  the  stator,  with  p^  ana 

P2  replacing  p2  and  p-j ,  respectively,  and  the  equation  is  solved  for  p  to 

l 

yield 


P3(0)  :  P2(0)  +  (p 2(9)  +  Pre?j  N 


2%(B) 


(11) 
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Radial  equilibrium  corrections  were  also  calculated  here  but  were  found  to 
be  negligible  because  the  exit  stator  flow  was  nearly  axial. 

This  equation  for  P^o)  was  evaluated  for  all  of  the  cases  considered 
above  and  the  results  are  shown  in  Figs.  3^-  through  36  for  the  180  deg  dis¬ 
tortion.  The  distribution  for  the  peak  pressure  condition  in  Fig.  3*+  may 
be  compared  directly  with  the  p2  distribution  shown  in  Fig.  15  for  the  same 
condition.  It  is  seen  that  the  overall  agreement  between  the  unsteady  pre¬ 


diction  and  the  measured  p^  distribution  is  not  as  good  as  the  agreement 
indicated  in  Fig.  15  for  the  p^  distribution  (in  the  regions  away  from  the 
screen  edges) .  A  detailed  examination  of  this  case  indicates  that  the  stator 


inlet  angle,  a^,  is  generally  less  than  50  deg  over  the  range  1O<0<1oG  deg. 

From  Fig.  33  this  yields  a  relatively  large  value  of  ^  and  from  Eq.  (11) 

s 

this  in  turn  produces  a  large  term  to  be  added  to  P2.  Over  this  same  0 


range  the  predicted  p^  overestimates  the  measured  p^  by  a  significant  margin 
while  the  predicted  p2  is,  on  average,  in  good  agreement  with  measured  values. 
Beyond  0  =  160  deg  ag  grows  larger,  reaching  a  peak  value  of  68  deg  in  the 
neighborhood  of  2OO<0<:26O  deg,  after  which  returns  to  a  value  of  approxi¬ 
mately  55  deg  behind  the  screen.  In  the  region  of  maximum  a 2  the  stator 
pressure  rise,  <1/  ,  becomes  very  small  and  the  additive  term  in  Eq.  (11) 
also  becomes  very  small.  In  this  same  region  the  predicted  p^  underestimates 
the  measured  p^  by  a  significant  margin  while  the  predicted  is  either  in 
ament  with  or  only  slightly  overestimates  the  measured  values.  It  can 


be  seen  in  Eq.  (38)  in  Appendix  I  that  is  directly  relatable  to  VA  through 
a  tangent  relationship.  Thus  becomes  large  in  the  undistorted  region, 
where  VA  is  large,  and  it  becomes  small  in  the  distorted  region,  where  VA 
is  small.  It  would  appear,  then,  that  the  increased  discrepancies  in  the 
stator  exit  static  pressures  may  be  traced  to  the  initial  assumption  that  the 
axial,  velocity  is  constant  through  the  stage  along  a  constant  9  path. 

In  Fig.  35  the  p^  distribution  for  the  midflow  condition  is  shown. 

Although  there  is  no  comparable  p^  distribution  shown  in  this  report  for 

this  case,  a  valid  comparison  may  be  made  with  the  0  distribution  in  Fig.  20. 

R 

This  time  the  calculation  of  p^  at  the  stator  exit  appears  to  produce  rela¬ 
tively  little  deterioration  in  the  agreement  between  theory  and  experiment 
relative  to  the  static  pressure  at  the  rotor  exit.  This  case  was  run  at  a 
midflaw  condition  which  is  removed  a  bit  from  the  high  stall  conditions. 
Consequently  the  VA  distribution  exhibits  smaller  deviations  from  the  mean 
and  the  variations  in  ot^  are  also  small.  Kence  the  additive  term  in 
Eq.  (38)  will  remain  relatively  constant. 

For  completeness,  Figs.  36,  37?  and  38  show  the  p^  distribution  for 
the  180  deg  distortion  near  stall,  the  135  deg  distortion  at  peak  pressure, 
and  the  9°  d.eg  distortion  near  stall.  In  general  the  agreement  between  the 
predicted  and  measured  pressures  at  the  stator  exit  for  all  five  cases  is 
comparable  to  the  agreement  between  the  predicted  and  measured  pressures  at 
the  rotor  exit,  with  minor  deviations  from  case  to  case.  However,  the 
results  shown  here  are  not  felt  to  be  satisfactory  in  their  present  form  to 
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represent  the  compressor  stage  and  should  not  be  used  for  predictions  of 
the  distorted  flow  response  of  subsequent  stages.  Note  that  the  remarks  made 
at  the  end  of  the  last  section  are  pertinent  here  also,  particularly  in  view 
of  the  dependence  of  p^  on  pg. 

Possible  Sources  of  Error 

There  are  many  possible  sources  of  error  in  this  work  which  could  not 
be  examined  during  the  first  year's  work.  Indeed,  mans’-  of  these  possible 
errors  could  not  be  anticipated  until  the  work  was  well  along.  Some  of  these 
are  being  scrutinized  as  part  of  our  current  investigation  and  have  con¬ 
tributed  to  the  modifications  made  in  the  fluid  dynamic  modeling  of  this 
problem.  This  will  be  outlined  briefly  in  the  next  section.  The  purpose 
of  this  section  is  to  list  these  possible  sources  of  error,  to  try  to  specu¬ 
late  on  what  their  impact  might  be  on  the  final  results,  and  to  indicate  what 
means  we  might  employ  to  eliminate  or  mitigate  the  effects  of  these  errors. 

Consider  first  the  list  of  assumptions  presented  in  the  first  section 
of  this  report.  At  the  inception  of  this  work  it  was  felt  that  these 
assumptions  were  necessary  to  permit  a  rational  beginning  to  be  made  in  this 
analysis.  This  judgment  is  borne  out  by  the  ability  of  the  analysis  to  yield 
the  agreement  that  it  does.  Now  that  a  tractable  formulation  exists  it  should 
be  possible  to  re-examine  each  assumption  and,  where  possible,  to  modify  the 
analysis  to  determine  the  sensitivity  of  the  results  to  changes  in  each 


assumption. 


Within  the  framework  of  the  experimental  data  used  for  comparison  herein. 
Assumption  1  concerning  the  thermodynamic  and.  inviscid  behavior  of  this  low- 
pressure  ratio  compressor  is  felt  to  be  valid.  Assumptions  2,  3>  and  5  imply  > 
in  part,  that  each  blade  passage  is  independent  of  every  other  blade  passage. 
This  may  have  contributed  to  the  sharp  excursions  observed  in  the  predicted 
pressure  rise  at  the  screen  edges  in  view  of  the  influence  that  each  blade 
exerts  on  the  neighboring  flow  field  and  the  fact  that  this  influence  is 
neglected  here.  Assumption  7  restricts  the  unsteady  effects  to  blade  aero¬ 
dynamic  response  only.  This  too  may  have  some  bearing  on  the  sudden  sharp 
excursions  mentioned  above.  Assumptions  4,  8,  and  9  are  interim  measures 
which  are  necessary  because  there  is  no  data  both  below  and  above  stall  for 
either  the  unsteady  lift  on  a  cascaded  blade  or  the  unsteady  drag  on  any 
blade.  Assumption  6  appears  to  be  reasonable  in  view  of  the  high  (0.7) 
hub- tip  ratio  of  the  annulus. 

In  virtually  all  cases  examined  in  this  report  a  considerable  portion 
of  the  circumference  operated  beyond  the  stall  point  and  within  this  region 
a  significant  portion  operated  beyond  the  limit  of  the  available  unsteady 
airfoil  data.  For  example,  for  the  180  deg  distortion  at  the  peak  pressure 
condition  it  can  be  seen  from  Fig.  11  that  75  percent  of  the  circumference 
was  above  stall,  and  over  35  percent  of  the  circumference  operated  beyond 
a  =  1.8  which  was  the  maximum  value  in  the  unsteady  data  tabulation.  It  can 
be  assumed  that  the  extremely  deep  penetration  of  the  rotor  blades  beyond 
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the  static  stall  angle  and  the  lack  of  unsteady  data  in  that  region  contri¬ 
buted  some  error  to  the  final  results.  At  the  other  extreme  it  was  necessary 
to  extrapolate  the  steady- state  rotor  characteristics  (see  Appendix  III)  to 
small  values  of  (large  J3-^'  for  some  portions  of  the  midflow  condition. 

This  too  could  have  contributed  to  our  error.  In  addition,  there  may  have 
been  insufficient  measuring  stations  in  the  experiment  in  the  regions  of 
steep  gradients  in  pressure  or  yaw  angle.  Within  the  scope  of  our  present 
contract  the  impact  of  these  possible  errors  can  only  be  guessed  at.  In 
each  case  the  only  rational  corrective  procedure  is  to  perform  experiments, 
either  into  the  unsteady  behavior  of  blades  and  cascades  in  deep  stall,  or 
into  the  steady  behavior  of  rotors  far  beyond  the  normal  operating  range. 

A  final  item  for  consideration  here  is  not,  strictly  speaking,  an 
error,  but  instead  represents  a  deficiency  in  the  theory  which  must  be 
overcome  in  the  future.  The  results  shown  earlier  use  measured  upstream 
flow  properties  obtained  at  the  inlet  plane  during  an  inlet  distortion  experi¬ 
ment.  These  upstream  conditions  are  strongly  dependent  on  the  presence  of 
the  blade  row,  as  shown  in  Fig.  39-  Here  we  have  plotted  the  measured  stall 
angle  parameter,  a,  versus  circumferential  position,  9,  for  the  peak  pres¬ 
sure  condition  with  a  180  deg  distortion.  The  dashed  line  in  Fig.  39 
represents  the  angle  of  attack  distribution  that  would  be  measured  immediately 
behind  the  distortion  screen  far  upstream  of  the  blade  row,  and  the  solid  line 
is  the  actual  angle  of  attack  distribution  at  the  blade  row  inlet  (repeated 
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from  Fig.  11) .  It  is  seen  that  the  blade  row  has  an  important  effect  on  the 
inlet  flow,  and  produces  a  further  distortion  of  large  magnitude  on  the  pre¬ 
viously  distorted  flow.  At  present  it  is  not  possible  to  predict  the  change 
in  the  distortion  pattern  as  it  comes  under  the  influence  of  the  blade  row, 
and  it  is  necessary  to  measure  the  upstream  conditions  during  an  actual  dis¬ 
tortion  experiment.  However,  the  ultimate  purpose  of  this  program  is  to 
provide  the  designer  with  the  necessary  tools  to  predict  compressor  stability 
due  to  inlet  distortion  without  recourse  to  test.  Hence  the  designer  should 
have  at  his  disposal  an  a  priori  method  to  specify  a  distortion  pattern,  pre¬ 
dict  its  behavior  as  it  approaches  the  blade  row,  and  predict  blade  row 
response  to  this  distortion.  Furthermore,  this  should  be  a  fully  coupled 
system  which  accounts  for  the  mutual  interactions  of  blade  row  response  to 
the  distortion  and  blade  row  influence  on  the  distortion.  Our  current 
efforts  are  being  directed  toward  exploring  this  objective  and  are  described 


briefly  below. 


CONCLUSIONS  AND  BRIEF  REVIEW  OF  FUTURE  WORK 


During  the  first  year's  work  reported  on  herein  it  has  been  shewn  that 
an  unsteady  model  based  on  unsteady  aerodynamic  data  from  an  isolated  airfoil 
could  lead  to  a  reasonable  representation  of  the  circumferential  pressure  dis¬ 
tribution  of  a  rotor  in  a  distorted  flow.  It  is  seen  that  in  general  the 
unsteady  prediction  is  in  better  agreement  with  experiment  than  the  quasi¬ 
steady  prediction,  but  that  there  are  regions  where  both  theories  are  in 
significant  disagreement  with  the  data.  Although  this  model  has  yielded 
encouraging  results,  it  is  recognized  that  it  depends  strongly  on  a  number 
of  intermediate  transformations  and  assumptions.  Furthermore,  we  have  con¬ 
cluded  that  it  is  difficult  to  improve  or  refine  the  results  using  physically 
rational  techniques  within  the  framework  of  the  existing  theory.  Hence,  we 
have  undertaken  to  reformulate  the  problem  in  more  fundamental  terms  involv¬ 
ing  the  direct  use  of  compressor  and  blade  row  parameters,  including  loss 
and  turning  characteristics.  This  procedure  is  described  briefly  below. 

A  simple  unsteady  model  based  on  cascade  loss  and  turning  correlations 
has  been  derived  to  predict  the  pressure  rise  across  a  rotor  operating  in  a 
distorted  inflow.  This  nonlinear  model  is  obtained  from  a  consideration  of 
the  equations  governing  the  passage  of  a  two-dimensional  incompressible  flow 
between  two  blades  of  a  compressor,  and  includes  provision  for  unsteadiness 
by  means  of  time  variations  in  both  passage  kinetic  energy  and  blade  row  loss. 

To  date  the  application  of  this  revised  theory  has  shown  significant  improvement 
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over  the  results  obtained  by  our  aerodynamic  force  procedure  described  above. 
In  addition  to  its  simplicity,  it  has  the  further  advantage  of  being  directly 
related  to  the  blade  row  through  the  blade  loss  and  turning  characteristics, 
rather  than  the  indirect  scaling  relationships  of  the  previous  aerodynamic 
model.  Accordingly,  this  model  is  being  adopted  for  all  subsequent  work 
being  done  on  this  contract. 

An  analysis  is  also  being  formulated  to  introduce  the  effects  of  the 
presence  of  the  blade  row  on  the  upstream  disturbance.  As  ultimately 
envisioned  this  represents  a  nonlinear,  time- dependent  analysis  of  the  duct 
flow  coupled  with  the  nonlinear,  time- dependent  blade  passage  analysis  which 
will  permit  us  to  solve  our  problem  for  large  disturbances  which  can  be 
steady  or  time  dependent.  This  is  a  formidable  problem  and  it  is  felt  that 
its  implementation  can  be  justified  only  by  demonstrating  the  tractability 
of  its  constituent  parts,  in  simplified  form.  The  first  part  consisting  of 
the  nonlinear  blade  passage  loss  analysis  has  already  been  implemented  and 
appears  to  give  excellent  results.  A  simplified  version  of  the  correspond¬ 
ing  nonlinear  flow  field  analysis  was  recently  attempted  and  is  described 
below.  A  numerical  procedure  was  devised  to  solve  the  time- dependent 
vcrticity  and  stream  function  equations,  subject  to  a  downstream  boundary 
condition  involving  the  pressure  change  across  a  blade  row.  The  problem 
was  simplified  by  assuming  the  blades  were  stationary  with  respect  to  the 


upstream  distortion  (i.e.,  the  cascade  problem)  and  by  eliminating  time 
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dependency  in  the  blade  passage  analysis.  Results  obtained  to  date  axe  very 
encouraging  and  appear  to  have  physically  realistic  behavior. 

Although  the  problem  described  in  this  section  is  simple  in  form  (i.e., 
a  solution  of  the  lightly  loaded  cascade  interaction  at  relatively  low  distor¬ 
tion  levels),  we  feel  it  contains  many  of  the  key  essential  elements  of  the 
complete  problem  which  should  be  solved  in  future  work.  First,  it  demon¬ 
strates  the  mutual  interaction  effects  of  the  upstream  distortion  vorticity 
and  the  bound  vorticity  of  the  blade  row.  Second,  although  the  parameters 
chosen  for  study  represented  modest  perturbations,  they  were  not  "small" 
perturbations  in  the  linear  sense,  and  the  solution  represented  a  nonlinear 
process.  Third,  from  an  operational  point  of  view,  the  numerical  solution 
was  stable,  and  we  feel  that  any  extension  of  this  solution  to  the  more 
complicated  case  of  a  moving  blade  row  can  be  maintained  stable.  Finally, 
it  was  a  time-  and  space- dependent  solution,  although  it  was  used  here  to 
solve  a  quasi- steady  problem.  Hence,  in  the  future  it  will  be  possible  to 
extend  this  procedure  further  to  fully  time- dependent  systems.  This  time- 
dependent  nonlinear  model  will  permit  us  to  specify  a  distortion  pattern  at 
the  inlet  of  a  compressor,  predict  its  behavior  as  it  approaches  the  compressor 
stage,  and  predict  compressor  stage  response  to  this  distortion. 
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% 
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FA 

Fa 
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H 
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Annular  area,  ft2,  Eq.  (42),  or  angular  velocity  parameter, 

Eq.  (64) 

Airfoil  eemicnurd,  ft,  Eq.  (93) 

Angular  acceleration  parameter,  Eq.  (65) 

Airfoil  chord,  ft,  Eq.  (10) 

Drag  coefficient,  Eq.  (55) 

Axial  force  coefficient,  Eq.  (10) 

Tangential  force  coefficient,  Eq.  (10) 

Lift  coefficient,  Eq.  (55) 

Normalized  lift  coefficient,  =  1  at  stall,  Eq.  (59) 

Drag,  lb,  Eq.  (55) 

Drag  differences,  Eqs.  (70),  (71) 

Axial  force,  lb,  Eq.  (7) 

Tangential  force,  lb,  Eq.  (7) 

Gravitational  acceleration,  =  32.16  ft/sec  ,  Eq.  (4l) 
Momentum,  lb- sec,  Eq.  I 17) 

Fully  stalled  lift  multiplier,  Eq.  (76) 

Blade  span,  ft,  Eq.  (12) 

Lift,  lb,  Eq.  (55) 

Lift  differences,  Eqs.  (72),  (73) 

Fluid  mass,  lb-sec  /ft,  Eq.  (12) 


LIST  OF  SYMBOLS  (Continued) 


Number  of  blades  or  blade  passages,  Eq.  (42) 

Wheel  speed,  rpm,  Eq.  (45) 

Corrected  wheel  speed  ratio,  Eq.  (45) 

Static  pressure,  lb/ft2  or  in.  l^o,  Eq.  (1) 

Total  pressure,  lb/ft2  or  in.  J^O,  Eq.  (l) 

Radius,  ft,  Eq.  (123) 

Time,  sec,  Eq.  (12) 

Temperature  ratio,  Eq.  (43) 

Wheel  speed,  ft/sec,  Eq.  (4) 

Relative  velocity,  ft/sec,  Eq.  (84) 

Axial  velocity,  ft/sec,  Eq.  (3) 

Absolute  velocity,  ft/sec,  Eq.  (2) 

Tangential  swirl  velocity,  ft/sec,  Eq.  (4) 

Annular  weight  flow,  lb/sec,  Eq.  (4l) 

Absolute  velocity  angle,  Eq.  (37),  or  angle  of  attack, 
Eq.  (6),  deg  or  rad 

Chord  angle,  deg  or  rad,  Eq.  (6) 

Yaw  angle,  deg  or  rad,  Eq.  (3) 

Relative  velocity  angle,  deg  or  rad,  Eq.  (5) 

Total  pressure  ratio,  Eq.  (44) 
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0 

Minimum  loss  or  minimum  drag 

p 

Potential  flow 

ref 

Atmospheric  reference 

R 

Rotor  blade 

!  s 

i 

Stall  or  stator 

!  ss 

Steady  state 

i 

1  u 

Unsteady 

1  « 

i 

Tangential 

; 

!  Superscripts 

(  ‘  ) 

First  time  derivative,  d/dt 

•  • 

(  ) 

2  i  2 

Second  time  derivative,  d  /dt 

(  )* 

First  circumferential  derivative,  d/d9 

(  )" 

: 

2  2 

Second  circumferential  derivative,  d  /d0 

i 
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The  flow  through  a  blade  row  has  been  the  subject  of  analysis  in  a  number 
of  basic  textbooks  (e.g.,  Refs.  4,  5,  6).  The  derivation  contained  herein  is 
simply  for  the  convenience  of  the  reader  and  will  serve  to  establish  a  con¬ 
sistent  set  of  notation  for  use  in  this  report. 

Consider  the  compressor  blade  row  shown  in  Fig.  3.  In  the  work  that 
follows,  each  equation  will  be  derived  for  a  quasi-steady  fluid  flow  through 
the  stage,  and  will  then  be  simplified  according  to  the  assumptions  governing 
the  present  analysis. 

First  consider  the  continuity  of  the  flow  through  the  rotor.  The  incre¬ 
mental  mass  entering  any  blade  passage  of  area  per  unit  time  is 


A  m. 
At 
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Ax, 
A  t 


(12) 


and  this  must  be  equal  to  the  mass  leaving  the  passage  at  the  exit  plane. 


:  PzT  *2Vfl 


where  x  is  measured  in  the  axial  direction  and  where  V  is  the  axial  velocity. 

n. 

If  the  passage  height  is  constant,  then  and  up°n  setting  these 

expressions  equal  we  have 


^tVA,  =  ^2va2  (14) 

If  the  further  assumption  of  incompressible  flow  with  negligible  temperature 
rise  through  the  blade  row  is  made,  then  Eq.  (14)  reduces  simply  to 

Pi-Pz'-P  to) 


(16) 


Note  that  although  VA  is  constant  axially  it  can  vary  circumferentially. 

Now  consider  the  axial  force  on  the  blade  row,  which  is  the  sum  of  the 
force  due  to  the  change  in  axial  momentum  across  the  blade  row  and  the 
axial  pressure  force  on  the  blade  row.  First,  the  axial  momentum  into  the 
blade  row  is 


and  out  of  the  blade  row  is 
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so  the  axial  force  due  to  momentum  change  is 
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for  constant  passage  height.  There  will  also  be  an  axial  force  due  to  the 
axial  pressure  gradient  acting  across  each  blade  passage,  or,  for ' constant 
passage  height, 


FAr=  ft  (p2“  p, ) 


The  total  axial  force  is  obtained  from  the  sum  of  Eqs.  (19)  and  (20), 


VtUp2-p,+  -  p  vA2  ) 


With  the  simplifications  of  Eqs.  (15)  and  (16)  this  reduces  to 


V  ^(p2-p,) 


Finally,  the  tangential  force  on  the  blade  row  will  be  derived.  The 
tangential  momentum  into  the  blade  row  is  given  by  the  product  of  the  incre¬ 
mental  mass  and  the  inlet  swirl  velocity  in  the  relative  frame  of  reference. 
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for  constant  passage  height. 

In  che  simplified  analysis  of  this  study  the  tangential  pressure  gradient 
will  he  zero,  as  explained  below.  If  a  control  volume  is  chosen  which  con¬ 
tains  exactly  one  blade  and  one  blade  passage,  then  in  this  simple,  quasi¬ 
steady  fluid  dynamic  model  the  pressure  distribution  on  one  streamwise  face 
of  the  volume  will  be  identical  to  the  pressure  distribution  on  the  other 
streamwise  face.  Hence,  the  tangential  pressure  gradient  is  zero  and  there 
is  no  pressure  force  in  the  tangential  direction.  Thus  Eq.  (25)  represents 
the  only  tangential  force  acting  on  the  system.  Again  using  the  simplifica¬ 
tions  of  Eqs.  (15)  and  (16),  Eq.  (25)  reduces  to 


r"’V\'v 


(26) 


An  immediate  objective  of  this  analysis  is  to  provide  the  necessary 
equations  to  predict  the  exit  static  and  total  pressure  distributions.  The 
static  pressure  is  immediately  obtained  from  Eq.  (22)  as 


1»1 


P  =6  +  -^4 

2  wl  r£ 


(27) 


To  obtain  the  total  pressure,  rewrite  the  steady  form  of  Bernoulli’s  equation 
(Eq.  (1))  at  the  rotor  exit  as 


Po?:P2+  v  py, 


In  this  equation  p?  is  already  known  from  Eq.  (27),  and  from  Eq.  (15)  p  will 
be  replaced  by  the  constant  density  p.  This  leaves  as  the  unknown 
quantity.  From  the  exit  velocity  triangle  in  Fig.  4, 


♦,u-y 


(29) 


with  =  VA  (from  Eq.  (16)),  and  given  by  Eq.  (26)  as 


V  :  V - -ii _ 

&Z  ft,  t  Up  vA 


It  will  be  convenient  to  deal  with  the  various  aerodynamic  forces  in 
coefficient  form.  Hence  we  define  these  coefficients  in  terms  of  the  inlet 
relative  dynamic  pressure  as 


h2 


Fs  =  ^-?V,2cfC 


Fn=  2  f7?  c'cFa 


With  these  substitutions  Eqs.  (27)  and.  (31)  become 


„  -_p  +  j^ 

2  1  2  r/c 


(32) 


(33) 


P02*  P2+  2  P 


v2+/u_v  +  _££^\2 

A  \  2(t /c)Va/ 


(34) 


To  complete  the  analysis  of  the  flow  through  the  stage,  we  now  take 
the  rotor  exit  variables  obtained  above  and  regard  them  as  inlet  variables 
for  the  stator.  The  analysis  is  simpler  here  because  the  nonmoving  stator 
operates  in  a  fixed  reference  frame.  Hence  the  intermediate  unsteady  blade 
force  analysis  is  not  needed,  and  a  simple  set  of  expressions  relating 
stator  pressure  rise  to  inlet  flow  is  all  that  is  necessary. 

As  shown  in  Fig.  33 »  the  steady- state  stator  characteristics  are  given 
herein  as  the  static  pressure  coefficient  for  the  stator. 


Vp,ef>N 


IT2 


(35) 


versus  the  stator  inlet  angle  ojg.  Hence,  to  solve  for  the  stage  exit  static 
pressure, 


U3 


'  Vf’2+('Vp-ef)5!+s 


(36) 


it  is  only  necessary  to  know  the  rotor  exit  static  pressure,  Pg,  and  the 
stator  inlet  angle,  a„.  The  latter  may  be  used  to  extract  ih_  data  from 

c  S 

Fig.  33  whereupon  Eq.  (36)  is  completely  determined. 

From  the  bottom  portion  of  Fig.  4, a 0  is  seen  to  be 


Mon''  ) 

l  u'vn2 ) 


where  VA  is  used  in  place  of  in  accordance  with  Eq.  (16).  When  Eq.  (30) 
is  substituted  for  and  Eq.  (32)  is  used,  the  formula  for  becomes 


a?  :  tan" 


v;J  cFo 
u-vn  +-'■■ -■& 

2  (f/c )  w. 


All  of  the  parameters  in  Eq.  (38)  have  been  previously  defined  and  are 


either  rotor  inlet  quantities  or  easily  calculated  rotor  blade  quantities. 
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(Ref.  7)  are  used  here, 


(1*3) 


y  Ttot' 

R  518.7 


*  =-Poi+Prtf  (44) 

2116 


where  TTQT  is  inlet  total  temperature  and  518.7  is  standard  temperature,  both 
in  °R,  and  where  PQ1  is  relative  inlet  total  pressure,  Pref  is  local  atmo¬ 
spheric  total  pressure,  and  2116  is  standard  atmospheric  pressure,  all  in 
lb/ft  .  The  corrected  wheel  speed  ratio,  N  is  given  by 


—  n/ y'tR 

N  = 


(N4/?R  'OES 


(45) 


where  N  is  actual  wheel  speed  and  the  subscript  Des  denotes  the  design  con¬ 
dition.  In  the  static  pressure  coefficient  equation  P2-P1  is  the  static 
pressure  rise  through  the  rotor,  p  is  the  relative  inlet  static  pressure,  and 
Pref  and  N  have  already  been  defined  above.  Note  that  in  these  equations  the 
terry erature  ratio  is  denoted  by  rather  than  by  9  to  avoid  confusion  with 
the  circumferential  position  around  the  compressor. 

In  the  experimental,  work  performed  at  UARL  for  the  Pratt  &  Whitney 
Division  of  United  Aircraft  Corporation,  which  was  described  briefly  above, 
the  inlet  total  temperature  ratio  was  approximately  one.  Hence,  upon  sub¬ 
stituting  Eqs.  (4l),  (42),  (44),  and  the  numerical  values  T^si,  g  =  32.16, 
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Eq.  (39)  becomes 


,  _  68,051  pVA 

(P0l+P,ef)N 


(46) 


2 

which  has  the  units  of  weight  flow  per  unit  area,  (lb/sec)/ft  .  The  measured 
compressor  characteristics  for  the  first  rotor  (and  first  stator)  of  the  UARL 
three- stage  test  compressor  is  shown  in  Fig.  7.  These  and  other  measured 
steady- state  quantities  for  the  first  rotor  are  listed  in  Table  I. 

To  compute  the  steady- state  forces  acting  on  the  rotor  blades,  solve 
Eq.  (40)  for  p2~P1  and  substitute  this  into  Eq.  (8)  to  obtain 


r  -  (p<  *  pref)NVRr/c 

T.  1  w  2 


±/>V' 


(47) 


which  is  the  axial  force  coefficient.  For  a  given  geometry  we  know  r/c,  for 

an  atmospheric  inlet  we  know  p,  and  from  Table  I  we  know  N  and  ^  .  The 

R 

remaining  quantities  are  easily  obtained  in  terms  of  the  given  characteristics. 
From  Eq.  (l),  for  the  absolute  frame  of  reference,  write 


Pi*  PT£t-  .  ik 


*01*  pref.  _  /  vqi  \ 

-JrPV,2  \  V,  / 


(48) 


and  from  the  inlet  velocity  triangle  of  Fig.  4, 


sin  p.  - 


VA 


(49) 


so  Eq.  (48)  becomes 


and  Eq.  (47)  can  be  rewritten  in  terms  of  known  or  given  steady- state 
quantities. 


The  steady- state  tangential  force  coefficient  is  most  easily  obtained  from 
a  combination  of  Eqs.  (26)  and  (32), 


2 — 
c  w 


(52) 


from  which  zhe  tangential  force  coefficient  can  be  calculated  in  terms  of 
known  or  given  steady- state  quantities. 


f i^atssBss^ssw^y^wCw.  *it.\?rvz:r- 


Finally,  if  we  define  the  lift  and  drag  coefficients  in  terms  of  inlet 
quantities, 


(55) 


then  use  of  Eqs.  (10)  and  the  inverted  form  of  Eqs.  (7)  yield 

CL  :  C^sin  /3|  +  Cp^cos  /3| 


Use  of  the  tabulated  steady- state  rotor  parameters  in  Eqs.  (51),  (5*0,  and 

(5 6)  yields  the  rotor  lift  and  drag  variation  with  inlet  angle,  0^.  These 

are  plotted  in  Fig.  40  versus  the  more  convenient  rotor  blade  angle  of 

attack,  a  ,  obtained  from  Eq.  (6).  Also  in  this  figure  a  subscript  R  has 

been  introduced  into  the  and  notation  to  denote  rotor  blade  quantities. 

Because  of  the  limited  range  of  the  data  available  the  a  range  is  restricted 

R 

to  values  greater  than  8  deg,  and  it  was  necessary  tc  extrapolate  the  lift 

linearly  towards  a  =  0.  Also,  the  calculated  drag  appears  to  go  to  zero 
R 

near  8  deg,  and  it  is  impossible  to  extrapolate  this  characteristic  to  any 
lower  angle  of  attack.  It  will  be  shown  later  that  this  behavior  can  be 
partially  compensated  for  by  our  scaling  procedure. 


A 
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Over  the  linear  range  of  each  curve  (i.e.,  for  o  — 1)  the  angle  of 
attack  at  any  point  can  be  replaced  by  the  corresponding  lift  coefficient, 
so  we  can  write 


°R  -- 


Clr  *  Cl 


RO 


'lrs  "  clro 


■=  C, 


( o-B  <  | ) 


(59) 


°~a  - 


(crA<  i ) 


(60) 


*  * 

wnere  Cj^  and  CL  are  the  rotor  and  airfoil  lift  ratios,  respectively, 
which  equal  unity  at  the  stall  point.  To  cause  the  two  lift  curves  in 
Fig.  4l  to  match  over  the  linear  range,  set 


aR  =  aA=  a 


(61) 


or 


Ct-R  ~Clro 
r  _  r 

“lrs  ^lro 


(62) 


A  comparison  between  Cj^  and  CL^  is  shown  in  the  upper  half  of  Fig.  42 

where  both  normalized  lift  curves  are  plotted  versus  a.  From  Fig.  40  the 

numerical  values  for  the  rotor  blade  are:  a  =12  deg,  CT  =  0.432: 

Rs  -tRs 

aRo  =  3*29  deg  (from  the  minimum  loss  condition  for  this  blade  row), 

CLro  =  0.265.  The  comparable  data  for  the  NACA  0012  isolated  airfoil  used 
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here  (Ref.  3)  are  ar^g  =  10  deg,  =  1.02.  The  definition  of  the  func- 

AS 

-ft 

tions  Cj^  and  Cj^  in  Eq.  (62)  forces  agreement  up  to  a  =  1,  and  it  is  seen 

in  Fig.  42  that  the  two  functions  continue  to  agree  up  to  a  =1.3*  Beyond 

this  point  the  two  curves  diverge.  To  obtain  congruence  over  the  entire 

range  of  o  it  is  necessary  to  introduce  a  correction  increment  6  in  the 

L 


Cl-t  =  CU  +  SL 


-ft  ft 

where  6L  is  defined  by  Eq.  (63)  as  the  difference  between  C-^  and  C^,  and 
is  shown  in  the  lower  half  of  Fig.  42. 

Consider  next  the  procedure  for  determining  the  unsteady  lift.  Use 
has  been  made  here  of  the  extensive  tabulation  of  unsteady  lift  obtained  at 
the  United  Aircraft  Research  Laboratories  (UARL)  on  an  NACA  0012  airfoil 
section.  These  data  are  fully  documented  in  Refs.  3  and  8,  and  the  unsteady 
lift  coefficient  is  tabulated  in  Ref.  8  as  a  function  of  the  angular  posi¬ 
tion,  or A,  the  angular  velocity  parameter, 


and  the  angular  acceleration  parameter. 


The  airfoil  angle  of  attack,  can  he  written  in  terms  of  the  rotor  blade 

angle  of  attack,  a  ,  by  using  Eqs.  (57),  (58),  and  (6l), 

K 


aA(6)  -  aAS  <t(8)  = 


aAS  (aR(&)  -  aR0) 


This  is  a  function  of  circumferential  position,  9  ,  so  the  time  derivatives 
become 


daA  dd 
Qfl'  69  ‘  dt 


flaA'  (9) 


(67) 


where  £2  =  do/dt  is  the  rotor  angular  velocity  and  the  prime  denotes 
differentiation  with  respect  to  G.  Once  these  derivatives  are  known  the 
unsteady  airfoil  lift  can  be  obtained  from  the  tabulation  of  Ref.  8.  A 
key  assumption  made  at  this  point  is  that  the  unsteady  rotor  blade  lift  is 
inci ementally  related  to  steady  rotor  blade  lift  in  the  same  manner  that 
the  unsteady  and  steady  airfoil  lift  are  related.  In  other  words,  we  assume 
that  Eq.  (63)  is  also  valid  for  unsteady  lift  response. 

It  is  convenient  at  this  point  to  provide  a  brief  synopsis  of  the  steps 
followed-  in  calculating  unsteady  lift.  A  similar  synopsis  will  be  given  in 
the  next  Appendix  for  unsteady  drag,  and  a  detailed  summary  of  the  computa¬ 
tion  procedure  will  be  presented  in  Appendix  V.  First  compute  the 
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circumferential  inlet  angle,  /^(e),  from  Eq.  (5),  using  the  various  inlet 

quantities  and  Bernoulli's  equation  (1).  Next  compute  rotor  blade  angle 

of  attack,  crR,  from  Eq.  (6),  and  the  stall  angle  ratio,  a  ,  from  Eq.  (57)* 

The  airfoil  angle  of  attack  is  given  by  Eq.  (66),  and  after  computing  the 

derivatives  in  Eqs.  (67)  and  (68),  determine  A  and  B  from  Eqs.  (64)  and 

(65).  Next  enter  the  tabulation  in  Ref.  8  and  interpolate  a  value  of  CL 

A 

which  is  then  normalized  to  Cj^  by  the  right-hand  portion  of  Eq.  (60) . 
Apply  the  incremented,  correction  of  Eq.  (63)  to  obtain  Cj^  and  then  solve 
the  right-heind  portion  of  Eq.  (59)  for  unsteady  rotor  lift  coefficient, 

Cj^.  Insert  this  veilue  together  with  the  unsteady  rotor  blade  drag 
coefficient  (see  Appendix  IV)  into  Eq.  (56)  and  invert  to  obtain  the  axial 


and  tangential  force  coefficients,  CF^  and  Cp^.  Finally,  substitute  these 
force  coefficients  into  Eqs.  (8)  and  (9)  to  compute  the  exit  static  and 


total  pressures  at  each  circumferential  station, 0  . 
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UNSTEADY  DRAG  HYPOTHESIS  AND  SCALING  PROCEDURES 

The  initial  procedures  leading  to  the  determination  of  the  steady- state 
rotor  drag  shown  in  Fig.  40  have  already  been  discussed  above.  In  this 
Appendix  we  will  describe  the  subsequent  steps  necessary  to  relate  steady- 
state  airfoil  and  rotor  blade  drag  coefficients,  and  further,  to  obtain 
unsteady  drag  from  steady- state  drag.  The  steady- state  drag  transformation 
procedure  is  similar  to  that  used  in  transforming  steady- state  lift. 

However,  in  view  of  the  complete  lack  of  unsteady  drag  information  it  was 
necessary  to  construct  a  hypothetical  unsteady  drag  model.  These  procedures 
are  described  below. 

.  irst  consider  the  steady-state  drag.  No  drag  data,  either  steady  or 
unsteady,  were  obtained  during  the  test  described  in  Ref.  8.  Hence,  steady 
data  (unpublished)  from  other  tests  conducted  at  UARL  on  an  NACA  0012  air¬ 
foil  will  be  used.  This  has  been  reduced  to  a  dependency  on  a  and  is 
presented  in  the  upper  half  of  Fig,  43.  Also  shown  for  comparison  is  the 
steady  rotor  blade  drag  from  Fig.  40,  also  plotted  as  a  function  of  a . 

As  before,  the  two  curves  v.ixl  be  made  congruent  by  computing  an  incre¬ 
mental  correction,  6^,  defined  by  the  equation 

CD.*C0.+  SD  i&9) 

n  A 
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where  6^  is  plotted  in  the  lower  half  of  Fig.  43. 

A  conceptual  model  will  now  he  constructed  to  permit  the  determination 
of  the  unsteady  drag.  This  model  is  based  on  two  well-known  facts;  namely, 
that  the  drag  rises  rapidly  when  the  flow  over  the  airfoil  separates,  and 
that  separation  is  delayed  on  an  airfoil  executing  an  unsteady  motion 
(Refs.  9,  10,  11).  The  technique  to  be  enployed  here  was  conceived  at 
UARL  and  is  described  in  Ref.  12.  The  basis  of  this  technique  is  the 
departure  of  the  instantaneous  lift  from  some  reference  value  of  the  lift. 

For  simplicity,  consider  first  the  steady- state  behavior,  as  shown  in 
the  sketch  in  Fig.  44.  In  the  lower  half  of  this  figure  the  solid  line 
represents  the  steady  drag  coefficient,  C^.  The  broken  line  labeled 
CDNS  is  hypothetical  extension  of  the  "no-stall"  drag;  i.e.,  we  presume 
that  if  no  stall  occurred  the  drag  would  continue  at  a  low  level,  as  extra¬ 
polated  here.  The  broken  line  labeled  Cnss  sin  a  represents  the  variation 
of  the  streamwise  component  of  the  normal  force.  At  the  point  where 
Cnss  sin  a  first  equals  the  steady- state  drag  we  assume  that  the  airfoil 
is  fully  stalled  and  we  label  the  angle  of  attack  accordingly  as  apg*  In 
the  upper  half  of  this  figure  the  solid  line  represents  the  steady-state 
lift  variation  with  a  ,  and  the  broken  line  labeled  Cr  represents  the 
straight  line  extrapolation  of  the  potential  flew  lift.  Suppose  we  have 
complete  lift  information  (i.e.,  the  entire  upper  half  of  the  figure)  but 
our  knowledge  of  the  drag  is  limited  to  the  no  stall  drag  curve,  C^g,  the 
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value  of  the  drag  at  a  fully  stalled  condition,  C^g,  and  the  angle  of 
attack  at  full  stall, «pS.  Furthermore,  suppose  we  wish  to  compute  the 
drag  variation  at  any  arbitrary  angle  of  attack,  a  .  The  following  hypo¬ 
thesis  is  proposed.  Let  be  the  unknown  departure  of  the  drag  from  the  no 
stall  drag,  and  let  Dg  be  the  known  difference  between  the  drag  at  full 
stall  and  “he  no  stall  drag  at  apg,  or 

Di  --  CDss  (o)  -  Cdns  (a)  (70) 


°2  =  C°FS  ’  C°NS  (aFs) 


Also  let  Lj_  be  the  known  departure  of  the  lift  from  the  potential  lift  at 
a  ,  and  let  Iv,  be  the  known  departure  of  the  lift  from  the  potential  lift 
at  a  fully  stalled  condition,  or 


L,  --  CLp  (a)  -  CL$S  (a) 


(72) 


CLp  (aFS  *  '  CL| 


Our  hypothesis  is  that  the  ratio  of  the  two  drag  departures  from  no  stall 
is  equal  to  the  ratio  of  the  two  lift  departures  from  potential  flow,  or 


IP! - " 


m 


When  we  substitute  Eq.  (70)  for  and  solve  for  the  unknown  steady-state 
drag  coefficient  we  obtain 

L,  (a) 

coss( a)  =  coNS(a)  +  ~ ~  °2  (75) 


This  hypothesis  was  tested  in  Ref.  12  and  the  results,  reproduced  from  Fig.  8 
of  this  reference,  are  shown  in  Fig.  45.  It  is  seen  that  the  agreement 
between  predicted  and  actual  CD  is  excellent. 

This  ability  of  the  hypothetical  drag  model  to  predict  steady- state 
drag  was  the  basis  for  the  unsteady  drag  hypothesis  which  follows.  Here 
we  assume  that  we  have  complete  steady- state  information  for  both  lift 
and  drag,  and  further,  that  we  are  given  an  unsteady  lift  hysteresis  loop. 

For  illustrative  purposes  we  shall  use  an  unsteady  lift  hysteresis  loop 
adapted  from  Ref.  3.  This  is  shown  at  the  solid  curve  denoted  by  CL  in 
the  upper  half  of  Fig.  46.  The  arrows  indicate  the  direction  of  increasing 
time.  The  dashed  line  is  the  steady  lift,  ,  the  dash-dot  line  is  the 
potential  flow  lift  curve  extension,  Cpp,  and  the  dotted  line  represents 
a  lower  limit  to  the  lift  characteristic  of  a  fully  stalled  condition,  CpFg. 
This  lower  limit  should  be  the  lower  envelope  (or  the  greatest  lower  bound) 
of  all  possible  lift  hysteresis  loops.  Although  it  cannot  be  determined 
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precisely,  experience  has  shown  that  this  bound  usually  occurs  at  approxi¬ 
mately  50  percent  of  the  steady  lift  at  any  angle  of  attack.  Hence  we  will 
define 

clfs  5  K  CLSS  (76) 

where  K  =  0.5  in  the  present  report.  The  corresponding  drag  curves  are 

found  in  the  lower  half  of  the  figure.  The  dashed  curve  is  the  steady  drag, 

Cdssj  'khe  dash- dot  curve  is  the  extrapolated  no- stall  drag,  Cp^,  and  the 

dotted  curve  is  the  equivalent  drag  of  a  fully  stalled  airfoil,  CL  sin  a. 

SS 

The  solid  curve  for  the  unsteady  drag  is  not  known  a  priori  and  its  deter¬ 
mination  is  the  object  in  the  present  analysis. 

The  unsteady  drag  hypothesis  to  be  made  here  is  analogous  to  the  steady 
drag  hypothesis  in  that  the  unknown  drag  deviation  from  a  given  datum  will 
be  assumed  to  be  proportional  to  the  known  lift  deviation  from  a  comparable 
datum.  To  establish  this  hypothesis  we  assume  that  (l)  if  the  lift  followed 
the  potential  flow  Cpp  curve  there  would  be  little  or  no  separation  and 
we  would  expect  that  the  drag  would  follow  the  no  stall  Cp^g  curve,  and 
(2)  if  the  lift  followed  the  fully  stalled  Cp^g  curve  this  would  constitute 
a  fully  separated  condition  and  the  drag  would  follow  its  fully  stalled 
counterpart,  Cpgg  sin  a.  Thus,  if  the  unsteady  lift  Cp^  is  between  the 
steady  and  the  potential  flow  lift,  we  assume  the  unsteady  drag  to  be  located 
proportionately  below  Cpgg  and  between  Cp  and  Cp^g,  or 
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C°u  =  c°ss 


.  /  Clu  ~  Clss  \ 
;  \  Clp  "  Clss/ 


Dee  '  CDNS) 


°Lu  “  Clss 
Clp -Clss  I  >« 


where  €  is  some  tolerance  limit.  Conversely,  if  the  unsteady  lift  Cj^  is 
between  the  steady  and  the  fully  stalled  lift,  we  assume  the  unsteady  drag 
to  be  located  proportionately  above  Cj^  and  between  Cp  and  Cj^  sin  a. 


:DU  =  CDSS  + 


/  CLSS"C|-U  \ 
\  ClSS  “  CLFS  / 


(Cm  sin  a  -  CD,J 


c«-u  "  C>-ss 

The  tolerance  restriction  in  Eq.  (77)  is  necessary  to  avoid  a  singular 
behavior.  If  in  fact  Cj^  is  nearly  equal  to  C^,  which  occurs  under 
unstalled  conditions,  we  will  usually  have  a  very  small  value  of  the  drag 
to  begin  with,  so  a  reasonable  replacement  for  Eq.  (77)  under  such  circum¬ 


stances  is 


Cqu  =  cDns "  cDns 


/ CLU  CLp  \ 


Clu  Clss 
|clp-clss|s  * 

In  the  present  analysis  the  tolerance  was  taken  to  be  €  =  10  . 

These  equations  were  used  to  calculate  the  unsteady  drag  coefficient 
shown  by  the  solid  curve  in  the  lower  half  of  Fig.  46.  As  before  the 


arrows  indicate  the  direction  of  increasing  time.  It  can  be  seen  that 
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the  hypothetical  unsteady  drag  loop  is  physically  reasonable.  In  regions 
of  enhanced  lift,  where  separation  is  delayed,  Cj^  is  smaller  than  the 
steady  drag,  and  in  regions  of  lift  deterioration,  where  the  flow  is 
heavily  separated,  Cj^  is  larger  than  the  steady  drag. 

This  procedure  will  now  be  summarized  briefly.  As  part  of  the  lift 
determination,  o  is  calculated  at  each  circumferential  position,  0 .  Use 
this  value  to  enter  Fig.  43  to  obtain  values  of  the  steady  airfoil  drag, 

Cda  =  and  the  drag  increment,  6^.  From  the  computed  unsteady  airfoil 

lift  and  the  steady- state  airfoil  characteristics  (previously  given)  compute 
unsteady  airfoil  drag  from  Eqs.  (77),  (78),  or  (79),  whichever  is  appro¬ 
priate.  Add  oD  to  this  value,  as  in  Eq.  (69),  to  obtain  unsteady  rotor 
blade  drag,  substitute  this  into  Eq.  (56)  together  with  the  unsteady  lift 
coefficient  obtained  previously,  invert  Eq.  (56),  end  proceed  with  the  cal¬ 


culation,  as  described  earlier. 


APPENDIX  V 

CALCULATION  PROCEDURES 

The  circumferential  position,  0 ,  is  the  primary  independent  variable 
in  this  computation,  and  each  of  the  following  steps  is  carried  out  at  each 
given  value  of  0.  In  the  calculations  performed  herein  the  input  0  incre¬ 
ments  were  chosen  to  produce  a  reasonable  representation  of  the  input 
variables  by  using  a  high  density  mesh  in  the  neighborhood  of  the  screen 
edges  and  a  low  density  mesh  away  from  the  screen  edges.  The  following 
summary  indicates  the  procedures  and  the  equations  used  to  calculate  the 
results  obtained. 

1.  Given  PO1(0),  p1(0),  produce  graphically  a  smoothed  pressure  differ¬ 
ence  variation  PQ1(0)  -  p^G),  introduce  this,  Poi(0)>  and  a  (0)  as  input. 

y 

2.  Confute 

v0,(9l  =Vf-"V9>-  P,<«»  (80) 

VA(©i  =  V0,ie)  sin  ayC8)  (8l) 

U  +VO|(0)co$ay(0)  (82) 
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...... 


t  I 

fc  , 


0{(B)  =  tan"'[vAl0)/Vfll0)] 


(83) 


V,(fl)  W  Va  (9)+  vo  (0) 


(84) 


aR(0)  =  /3* ~  /3,(0) 


(85) 


where 


where 


£*«  46.29  deg 


aD(0)-  ac 


<Vfl)  =  a„  -a 


RS  Ro 


a  s  3. 29  deg 


a  =  12  deg 
RS 


o(0)  =  a  a(0) 
A  As 


(86) 


(87) 


(88) 


(89) 


where 


(90) 


a  z  10  deg 


3.  Use  a  numerical  differentiation  procedure  (such  as  that  described 
in  Ref.  13)  to  compute 


a'(6)  =  da.  /  d0 
A  A 


(9D 


a"  (9)  =  d2aA/d29 


(92) 


and  from  the  given  values  of  b  and  ,  compute 


a(0)  = 


V,  (0) 


(93) 


□  (0)-  b^2/A(0) 
B(0)'  \f(0T 


(9*0 


b.  At  each  value  of  0  use  a  ,  A,  and  B  to  enter  the  existing  tab,.la- 

A 


tion  of  unsteady  lift  and  interpolate  to  obtain  the  local  value  of  the 

unsteady  airfoil  lift,  CT  ,  and  compute  the  unsteady  normalized  lift  from 

Mi 


C,  ( 9 )  =  C.  (0)  /  1.02 


Lu 


(95) 


Cb 


(as  in  Eq.  (62)  for  steady-state  lift).  At  this  same  value  of  a  ,  obtain 

Pl 

the  steady-state  airfoil  lift  by  interpolating  at  A  =  B  =  0,  and  call  it 


C 


L 


(6)  =  c,  (0) 

ss  la 


(A=B=0) 


(96) 


Also  compute 


C.  (0 )  =  K C  (0) 

lFS  lSS 


where  K  will  probably  be  K  =  0.5,  and  compute 


C,  (0)  sina  (0) 

Lcc  A 


(97) 


(98) 


5.  Next  we  compute  the  airfoil  drag  coefficient.  Enter  Fig.  43  at 
the  appropriate  value  of  o  for  each  n  and  read  the  airfoil  steady- state  drag, 
calling  it 


C 


D 


(0)  =  Cn(0) 

SS  UA 


(99) 


Test  Cj^  versus  C^gg,  and  check  the  difference  between  Cjgp  and  Cj^g.  If 

C,  (0)>C,  (0)  (100' 

Lu  Lss 

and  if 
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aTO!R5^rr=^a^?;  rzz7&z??m W 


tc,  iewc.  teu  >« 

Lp  uss 


(101) 


(where  €  is  some  tolerance, 


taken  to  be  e  =  10" 2  herein)  then  compute 


cje)-cnje) 


_/Ct-ge)~Ct-s,^\(c  (g)- 

\C40)~C*4O)/  ss 


C0  (0)) 

UNS 


(102) 


Here  and  in  the 


next  equation  for  C„  .  the  no-stall  <Jro?  ooeffioient,  has 


been  idealized  as  two  straight  line  segments. 


.0025cr+  .008,  0<tr<  8 
. 0055556 (cr-H),  .8^  ^  -  2.6 


(103) 


C.  (0)*  C.  {6) 

Lu  Lss 


(104) 


and  if 


1 C,  (61-  C,  (0)1  S  e 


(105) 


then  compute 


„  /C,(0)-C  ,J9)\ 

cn<e)  =  cD  (0)-c  <01  I  Lu  L.g— 

0(1  dns  uns  \  C.  (0)  / 

\  UP  t 


(106) 
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Finally,  if 
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The  normalized  lift  in  Eqs.  (109)  and  (111)  are  converted  to  rotor  lift  by- 


solving  Eq.  (62)  for  Cj^(o)  and  inserting  the  appropriate  numerical  values 
of  the  constants.  In  the  present  case,  for  both  steady  and  unsteady  lift, 
the  formula  is 


C.  (0)  =  .167 C ,  (0)  +  .265 
lr  lr 


(113) 


7.  Next  compute  the  steady  and  the  unsteady  force  coefficients.  Use 
Eq.  (7)  written  in  coefficient  form. 


C_(fl)  =  C.  (0)  cos  /3A9)  -  Cft(0)sin/3  (0) 
''a  lr  ur 


(114) 


Cc(0)  =  C,  (0)sin/3.(0)  +  C.  (0)cos/3  (0) 
lr  dR 


(115) 


8.  Now  compute  the  static  and  total  pressure  profiles  at  the  rotor 
exit  plane,  using  Eqs.  (8)  and  (9), 


2 

V8'  *  V«>  +  -^7T  CF.18’ 


(116) 


PO2(0)  •  p2(0) 


I  f  2  /  V|2(0)CW01 


(117) 


9.  From  Eq.  (40)  compute  the  static  pressure  coefficient, 


p,(0)-pt(0) 


R  (p.Cei+P^N1 


(116) 


and  from  Eq.  (46)  compute  the  flow  coefficient 


<p(B)  = 


68,051  pVA(0) 
(P0,(8)+Pr,f)N 


(U9) 


n 

where  Pref  is  the  local  atmospheric  pressure  in  pounds/ft  .  The  corrected 
wheel  speed  ratio  is  defined  in  Eq.  (45)  and  was  numerically  given  by 


N  =  0-4495 


(-120) 


in  this  analysis. 

10.  To  calculate  stator  exit  quantities,  first  compute  the  stator 


inlet  angle.  '  From  Eq.  (38)  write 


a  =  tan 
2 


V®) 


v,  (0)CF  ( 9 ) 


(121) 


I  U— Vrt{0)  + - —  1 

V  2(t/c)Va(0)/ 

The  denominator  in  this  equation  was  previously  computed  as  part  of 
Eq.  (117).  At  each  value  of  0,  enter  Fig.  33  at  the  computed  value  of 
a  2  and  obtain  a  value  for  Then  use  Eq.  (11)  to  calculate  the  stage 

exit  static  pressure  variation, 


p3(0)  =  P 2<0)+  (p2(0)  +  Pr#f)N2*  (0) 


(122) 
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APPENDIX  VI 


RADIAL  EQUILIBRIUM  CORRECTION 


If  there  is  no  radial  component  of  velocity,  the  fluid  is  in  radial 
equilibrium  and  the  equation  governing  the  static  pressure  can  be  simplified 
to  (cf.  p.  8  of  Ref.  6) 


dp  2 

IT  *  pu>  T 


(123) 


where  o  is  the  (constant)  angular  velocity  of  the  fluid  in  the  absolute 
frame  of  reference.  When  Eq.  (123)  is  integrated  between  the  mid-annulus 
(subscript  M)  and  the  wall  (subscript  W)  the  result  is 


pw  pm  =  2  pu  lrw  rM) 


(124) 


In  Fig.  4  it  is  seen  that  the  absolute  tangential  velocity  at  the  exit  plane 
is  VpJ tan  a^,  and  at  the  mid-annulus  this  is  related  to  u  by  the  formula 


fe)M- 


(125) 


Thus,  the  mid-annulus  pressure  from  Eq.  (124)  is  reduced  from  the  wall  value 
according  to  the  equation 


pw  2  p(  ton\)M[(7^) 


(126) 
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TABLE  I 

STEADY- STATE  ROTOR  CHARACTERISTICS 
(Axial  Flow  Inlet,  a  ,  =  90  deg) 


lb/ sec 


ft 


12.8 

13.2 

13-7 

14.1 
14.55 

15.3 

16.4 

17.4 

18.1 

19.2 


0.0243 

0.0266 

0.0288 

0.0305 

0.0320 

0.0338 

O.0352 

0.0355 

0.0353 

0.0347 


a2 

VA 

V01 

v 

v02 

POI  + 
Pref 

! 

deg 

ft/sec 

ft/sec 

ft/sec 

lb/ft2 

43.9 

76.0 

76.1 

109.5 

2110.93 

45.6 

78.3 

78.4 

109.4 

2110.52 

47.8 

81.3 

81.4 

109.6 

2110.26 

49.5 

83.8 

83.9 

109.9 

2109.97 

51.5 

86.6 

86.7 

110.4 

2109.72 

54.5 

90.9 

91.0 

111.3 

2109.33 

59-0 

97.5 

97.7 

113.5 

2108.86 

63.O 

103.4 

103.6 

115.9 

2108.48 

65.6 

107.6 

107.8 

118.0 

2108.19 

69.6 

114.2 

114.4 

121.6 

2107.79 

0.4512 

0.4503 

0.4503 

0.4503 

0.4501 

0.4503 

0.4506 

0.4504 

0.4503 

0.4501 
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MEASURED  CIRCUMFERENTIAL  TOTAL  PRESSURE  PROFILE 
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-  PEAK  PRESSURE  CONDITION 
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